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Abstract 



A sample of 2.2 million hadronic Z decays, selected from the data recorded 
by the Delphi detector at Lep during 1994-1995 was used for an improved 
measurement of inclusive distributions of , and p and their antiparticles 
in gluon and quark jets. The production spectra of the individual identified 
particles were found to be softer in gluon jets compared to quark jets, with a 
higher multiplicity in gluon jets as observed for inclusive charged particles. A 
significant proton enhancement in gluon jets is observed indicating that baryon 
production proceeds directly from colour objects. The maxima, of the ^- 
distributions for kaons in gluon and quark jets are observed to be different. 
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1 Introduction 

The different colour charges of quarks and gluons lead to specific differences in the par- 
ticle multiplicity, the energy spectrum and the angular distributions of the corresponding 
jets. Beyond the study of these differences [Q, which are related to the perturbative 
properties of QCDQ elementary fields, the comparison of gluon and quark jets opens up 
the possibility to infer properties of the non-perturbative formation of hadrons directly. 

The study of ratios of identified {n"^, K^,p{p)) particle distributions in gluon (^f) and 
quark (g) jets is the main subject of this paper. 

Gluon jets are selected in bbg events by tagging the h quarks using techniques based 
on the large impact parameters of tracks coming from heavy particle decays. The Ring 
Imaging Cherenkov Counters (RICH) of the Delphi detector provide particle identifica- 
tion over a wide momentum range in combination with the ionization loss measurement 
of the Time Projection Chamber (TPC) and so allow a detailed comparison of identified 
particle spectra in gluon and quark jets. These are used for a detailed test of QCD based 
fragmentation models and also to check MLLA0 and LPIID[| predictions 0. 

This paper is organized as follows. In Section ^ the hadronic event selection, the 
quark/gluon separation, and the particle identification are described briefly. The experi- 
mental results are presented and compared with the predictions of models in Section ^. 
Finally a summary and conclusions are presented in Section ^. 

2 Experimental Technique and Event Sample 

A description of the Delphi detector, together with a description of its performance, 
can be found in 

2.1 Event Selections 

The data collected by Delphi during 1994-1995 are considered in the present analysis, 
during which time the RICH detectors (the main particle identification detectors) were 
fully operational and the Vertex detector was equipped with a three-dimensional readout. 
The cuts applied to charged and neutral particles and to events in order to select hadronic 
Z decays are identical to those given in ^ and [^]. The data sample passing the selection 
of hadronic events contained 1,775,230 events with a small contamination (< 0.7%) arising 
from T^T~ pairs, beam-gas scattering and 77 interactions |0. 

The influence of the detector performance on the analysis was studied with the full 
Delphi simulation program, Delsim [Q. Events generated with the Jetset 7.3 Parton 
Shower (PS) model [^, with parameters tuned by Delphi 0, were passed through 
Delsim and processed with the same reconstruction and analysis programs as the real 
data. 

Three-jet events were clustered using the Durham algorithm ||^ with a jet resolution 
parameter Ucut = 0.015. The value used for the cut-off was optimized using the Jetset 
7.3 PS model, by maximizing the statistics available and the quark/gluon purity attained 
for the three-jet event samples 0. 

The jet axes were projected onto the event plane, defined as the plane perpendicular 
to the smallest sphericity eigenvector obtained from the quadratic momentum tensor, 

^ QuantumChromoDynamics 
^Modified Leading Log Approximation 
^ Local Parton Hadron Duality 
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— Yl^=iPiaPii3- The jets Were numbered in decreasing order of jet energy, where 
the energy of each jet is calculated from the angles between the jets assuming massless 
kinematics: 

sin6'i + sin6'2 + sin6'3 ^ ' ' ' ( ) 

where 6j is the interjet angle as defined in Figure |1|. 

For a detailed comparison of quark and gluon jet properties, it is necessary to ob- 
tain samples of quark and gluon jets with similar kinematics and the same underlying 
scales |]IU| . To fulfill this condition, two different event topologies were used, as illustrated 



in Figure 0: 

mirror symmetric events, with 62 and 63 G [150° — 15°, 150° + 15°], 
subsequently called Y events, and 

three-fold symmetric events, with 62 and 63 G [120° - 15°, 120° + 15° 
subsequently called Mercedes events. 




Figure 1: Event topologies of symmetric Y events and Mercedes events; 
6i are the angles between the jets after projection into the event plane. 

For Y events only the low energy jets (jets 2 and 3 in Figure |lD were used in the analy- 
sis. For Mercedes events all jets were used in the analysis. The appropriate scale for these 
jets, equivalent to the e~^e~ beam energy can be approximated by k = Ejetsm9i/2 [1T0| . 
Mercedes events are mainly used to study the scale dependence of particle production. 

In order to enhance the contribution from events with three well-defined jets attributed 
to qqg production, further cuts (sum of angles between jet, polar angle of each jet axis, 
visible jet energy per jet and number of particles in each jet) are applied to the three-jet 
event samples, as in [^]. The number of three-jet events in the Mercedes and Y samples 
is 11,685 and 110,628 respectively. 



2.2 Quark and Gluon Jet Identification 

The identification of gluon jets by anti-tagging of heavy quark jets is identical to that 
described in [^. Heavy quark tagging is based on large impact parameters with respect 
to the primary vertex due to the long lifetime of the heavy particles. 
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The efficiency and purity calculations were made using events generated by the Jetset 
7.3 Monte Carlo model tuned to Delphi data and passed through Delsim. 

Even in simulated events, the assignment of parton flavours to the jets is not unique, 
as the decay history is interrupted by the building of strings in models such as Jetset 
or by the parton assignment of clusters in the case of Herwig. Thus two independent 
ways of defining the gluon jet in the fully simulated events were investigated. The first 
method assumed that the jet which has the largest angle to hadrons containing heavy 
quarks is the gluon induced jet Q (angle assignment) and in the second method the jet 
containing the fewest decay particles from the heavy hadrons was assigned to the gluon 
(history assignment). Both methods give similar results and therefore the purities can be 
estimated with small systematic uncertainties P]. 

Gluon jet purities of ~ 82% for Y events and Mercedes events were achieved. Here 
the purity is defined as the ratio of correctly identified gluon jets to the total number of 
jets tagged as gluons. There are 24,449 events with an identified gluon jet in the case of 
Y events and 1,806 in the case of Mercedes events. 

2.2.1 Corrections 



Y 


quark content (dusc) 


b quark content 


gluon content 


normal mixture 


49.5% 


1.6% 


48.9% 


b tagged jets 


25.1% 


58.2% 


16.6% 


gluon tagged jets 


13.7% 


4.2% 


82.0% 


Mercedes 


quark content (dusc) 


b quark content 


gluon content 


normal mixture 


64.2% 


2.3% 


33.3% 


b tagged jets 


17.3% 


73.6% 


9.1% 


gluon tagged jets 


11.1% 


6.8% 


81.8% 



Table 1: Compositions of different jet classes in Y and Mercedes events. The statistical 
errors are smaller than 1%. 



Table |1| shows the fractions of "light" quark, b quark and gluon jets in the three 
different jet classes entering the analysis of Y and Mercedes events. The classes are 
normal mixture jets, gluon tagged jets and b-tagged jets. "Light" quark denotes here a 
mixture of dus and c quarks. The jets of the normal mixture are taken from events in 
which the heavy hadron tag failed. Therefore they are predominantly unidentified dusc 
quark and gluon jets. Denoting a data bin of an observable of a pure gluon, light or b 
quark jet sample with Rg, Rdusc and Rb respectively, the measured observables in the 
three tagged classes Rgtag, Rhag and Rmix can be written as: 



Rmix 


dusc 

Pmix 


Rdusc + Pmix 


Rb 


~^ Pmix 


■ Rg 


Rbtag 


dusc 

~ Pbtag 


Rdusc + Pbtag 


Rb 


+ PLg 


Rg 


Rgtag 


dusc 

~ Pgtag 


Rdusc + Pgt^g 


Rb 


+ 

^ ^gtag 


■ Rg 



where the are the fractions as e.g. shown in Table 0. The observables for the pure 
samples can then be obtained by solving equation ^ for Rg, Ri and Rb- The statistical 

"^There are almost always only two heavy hadrons in an event, because the g —t qq splitting into heavy quarks is strongly 
suppressed. 
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errors on the fractions p'^ are less than 1% and are fully propagated with only a small effect 
on the total errors. Instead of the "light" quark sample containing dus and c quarks, a 
general quark sample containing dusc and h quarks can be deduced by setting pl^^^ = 

and adding p^j^. to p^fx Pgtag Pg^ag^ reducing equation 2 to a 2 x 2 matrix equation. 
In a similar way the light quark sample can be reduced to only containing dus but no 
c quarks by hardening the heavy hadron tag and treating the c quarks like b quarks in 
the Monte-Carlo. This leads to the pure observables Rg, Rdus and Rch which are related 
to the measured observables Rgtagi Rhtag ^i^d Rmix in the same way as before but with 
different pj. To correct for the limited detector acceptance, secondary reinteraction of 
particles and resolution of the detector, an acceptance correction factor 

jdMC 
J^MC+detector 

is also applied to the data bin by bin for each distribution. Here R'^'^ denotes pure model 
distributions (referring to "light" quark and gluon jets) and T^^c'+'^etector ^jg^otes the full 
simulation including detector effects treated like the data. Long lived particles like the 
and the A° were considered as instable when computing model distributions. 



2.3 Identification of Final State Particles 

For the measurement of the 7r+ , and proton content in jets a combined tagging 
procedure based on the Cherenkov angle measurement in the RICH detector and on the 
ionization energy loss (dE/dx) in the TPC was applied which is described in detail in 0. 

The combined application of TPC and RICH allows a continuous particle identification 
in the momentum range of 0.3-45.0 GeV/c. Table ^ shows which detectors were used to 
identify pions, kaons, and protons depending on their momentum. 





Momentum Range [GeV/c] 


0.3 - 0.7 


0.7-0.9 


0.9 - 1.3 


1.3 - 2.7 


2.7-9.0 9.0-16.0 


16.0 - 45.0 


7r 


TPC 


LRICH S 


GRICH S 


K 


TPC 


LRICH S 


GRICH V 

+ 

LRICH S 


GRICH S 


P 


TPC 


TPC 

+ 

LRICH V 


LRICH S 


GRICH V 

+ 

LRICH S 


GRICH V 


GRICH S 



TPC Identification by measurement of the energy loss 

LRICH S(V) Signal (Veto)-Identification with the liquid RICH 
GRICH S(V) Signal (Veto)-Identification with the gas RICH 



Table 2: Application ranges of the detectors for particle identification 



An algorithm was developed to obtain an optimal combination of the particle iden- 
tification possibilities of the TPC and the RICH. It combines the probabilities for the 
particle identification with the TPC and the RICH by a simple multiplication and renor- 
malization, and predefines three different identification classes, loose, standard, and tight. 
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by using well chosen cuts on this combined probability distribution. These cuts for the 
particle identification probabilities allow particle identification performances with differ- 
ent purities R and efficiencies e: 



# of particles of kind i identified as kind j 

^ of all particles identified as kind j 
^ of particles of kind i identified as kind j 
7^ of all particles of kind i 



Figure ^ shows the efficiency of the combined particle identification of pions, kaons, 
and protons as a function of the momentum of the particle. The curves of the expected 
energy loss and the Cherenkov angle, 6c, are shown in the upper part of Figure ^ 

Figure ^ shows the resulting purities of the particle identification for Y events. To keep 
the infiuence of particle reinteractions in the detector material small, this distribution is 
restricted to negatively charged particles in the momentum range p < 2.7GeV/c (for 
details see P]). The purity matrix is predominantly diagonal (Figure |^,f,k). The most 
important, however still negligible background for the pion reconstruction stems from 
electrons and muons. An exception are energetic electrons and muons from semilep- 
tonic hadron decays in b (and c) jets with an identification rate up to 20%. The main 
background for the kaon selection are pions. A kaon identification purity of ~ 70% is 
achieved. As the kaon production rate is almost one order of magnitude smaller than the 
pion production rate in hadronic Z decays, this implies a very efficient pion suppression. 
The lower proton identification purity in b jets is mainly due to the higher probability to 
identify kaons as protons because of the higher K multiplicity in b jets. 



Acceptance correction of the spectra of identified particles 



From the measured particle spectra In, Ik, and Ip of identified particles one obtains 
the spectra of pure hadrons S^, Sk, and 5"^ by solving the equation system: 



/ In 






( 




F^ 
^K 






Ik 








6^ 




^p 




\ Ip 


J 






^7r 


e^K 


^p 


1 



I Sn\ 

Sk 



(4) 



This correction is applied before the correction of the jet purity. The values e\ denote 
the efficiencies that the particles i are identified correctly; the values with i ^ j are 
proportional to the background of particle class j. A correction for secondary reinter- 
action of particles in the detector material was included in the overall correction factor 
Equation S. 



3 Results 

For identified particles in quark and gluon jets the multiplicty and the semi-inclusive 
distributions as function of the momentum p, = Inl/xp = Inpjet/p and the rapidity 

1] = In with respect to the jet axis have been measured. For each particle also the 

Ratio between gluon and quark jets R and the normalized ratio R' = R/tch of each of 
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p [GeV/c] 

Figure 2: Curves of expected values and efficiencies of the particle identification; a) shows 
the curve of the expected values for specific ionization for pions (vr), kaons {K), protons 
(p), muons (/i), and electrons (e) as a function of the momentum, b) shows the curve 
of the expected values for the Cherenkov angle 9c in the liquid and gas radiator for 
the same particle hypotheses. 6c'^°'^ was multiplied by a factor 5. The curves begin at 
p = 0.3GeV/c for the liquid radiator and at 1.7GeV/c for the gas radiator, c) shows 
the resulting efficiencies for Y events for the standard identification of pions, kaons and 
protons in the barrel of DELPHI for the 1994-95 data. Light vertical lines in all plots 
indicate the threshold-momenta of vr, K and p identification in the two RICH radiators. 
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1 1 ■ ■ ■ 


: i) 

: La.** * 


: j) 
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:<> 




: 1) 
: (1 

■ o 



1+ 



1+ 



^1 
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1 10 



10 



1 10 
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Figure 3: Purity of the particle identification in Y events. Here 'Gen.' denotes the generated flavour 
of the particle and 'Id.' denotes the tagged particle flavour, -i means not. 
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the observables have been studied, where r^h denotes the corresponding ratio obtained 
for all charged particles. 

Note that the particle multiplicity of jets is not a well defined subject which depends on 
details of the jet definition influencing the assignment of low momentum particles to the 
jets. The given results on multiplicities and also the particle distribution corresponding 
to very small momenta therefore always refer to the jet definition specified in Section |2.1| . 

Special emphasis here lies on the measurement of ratios in gluon to quark jets R as 
in these ratios the systematic error is considerably reduced as most of the systematic 
uncertainties cancel out. The double ratios R stress particle specific differences between 
gluon and quark jets. 

In this analysis gluon jets are in general compared to a "duscb" quark jet reference 
sample. The flavour mix of this sample is that of hadronic Z decays. For the comparison 
of particle multiplicities also reference samples were used where the b events ("dusc") 
and all heavy quark events (" dus" ) were removed. 



3.1 Multiplicities 

In Figure ^ we present the mean multiplicities Nq and Ng for identified particles in 
quark and gluon jets respectively, as well as their ratio, R = Ng/Nq, and the normalized 
multiphcity ratios, R' = R/tch- 
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• DELPHI, Quark=duscb 

O Jetset 7.4, def. bary. prod. 

▲ Jetset 7.4 

A Jetset 7.3 

■ Ariadne 4.06 

□ Herwig5.8C 



0.5 0.75 1 1.25 1.5 1.75 



^^'x- ^^X / ''ch 



Figure 4: Mean multiplicities A^^and ratios of multiplicities Rx for identified particles 
X in quark and gluon jets of Y events compared to different Monte Carlo models. 
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The multiplicities measured in quark jets for identified hadrons and for all charged 
hadrons depend on the composition of the quark flavours within the quark jet sample. 
The values obtained for the multiplicities A^^ and Nq and the ratios R resp. R' are given 
in Table ^. The determination of the systematic errors is described below. In Table ^ the 
normalized multiplicity ratios R' are compared to the predictions from the Monte Carlo 
simulations for Y and Mercedes events. The data show a significant proton enhancement 
in gluon jets for Y events. A similar enhancement, although less significant, is also seen 
in Mercedes events. The slight change observed for i?^ and R'^ (see Table ^ for different 
flavour compositions can be understood due to a stronger K production and a depleted 
proton production in events with heavy quarks. 

Simulations with statistics superior to the data based on the Jetset 7.4 PS model. 



Ariadne 4.08 ||ril, and Herwig 5.8 ||T2|| with parameters tuned by Delphi [g are 
compared to the data. Herwig underestimates both the kaon and the proton production 
in gluon jets. In contrast Jetset and Ariadne tend to overestimate the proton 
production in gluon jets. The Jetset model with default baryon production f\ deviates 
less. The difference to the other model is that here the extra suppression at the string 
end, which had been introduced to describe baryon production at large scaled momenta 
better 0, is inactive. The excess of baryon production in gluon jets indicates that baryons 
are directly produced from a colour string and not via intermediate colour and baryon 
number neutral clusters. This is discussed in more detail in Section |3!4. 

As a cross-check the summed multiplicity ratio Rtt±+p±+k± was calculated. A value of 
1.21 ± 0.01 was obtained in the case of Y events and 1.29 ± 0.02 in the case of Mercedes 
events. Both numbers are in good agreement with a direct measurement of this ratio |T^ 
(Y: 1.235 ± 0.030, Mercedes: 1.276 ± 0.059). 

For completeness Table |^ shows a comparison with measurements of other exper- 
iments. A significant excess of proton production was observed by Argus and 
Opal |T^. No quantitative comparison is, however, possible due to the different energies 
or event topologies. 



Systematic Errors 

Table ^ summarizes the influences of the most important sources of systematic error 
for the determination of the multiplicities and their ratios. To obtain systematic errors 
comparable with the statistical errors, half the difference of the value obtained when a 
parameter is modified from its central value is quoted as the systematic uncertainty. The 
single errors are added quadratically. The following sources of systematic uncertainties 
were examined. 

1. Decays of K°,A° 

It was examined whether the ratios of the production rates of pions, kaons and pro- 
tons in quark and gluon jets are influenced by Kg and A° decays. These decays are 
reconstructed with the program MAMMOTH [jl5l for the data and detector simula- 



tion. At the generated level of the simulation, these particles have been treated as 
stable particles. 
2. Secondary interactions 

Another source of uncertainty stems from particles produced in reinteractions of 
primary particles with the detector material. Positively charged pions and protons 

^ Note that Ariadne employs the non-perturbativc hadronization model of Jetset. 
''diflferent from mstj{12) = 2 
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Y events 


Kind 


Quark 


Nquark Nciuon 


Rx R'x 


TT 


duscb 
dusc 
dus 


5.852 ± 0.036 ± 0.071 7.067 ± 0.032 ± 0.077 
5.702 ± 0.039 ± 0.069 7.043 ± 0.031 ± 0.076 
5.672 ± 0.040 ± 0.067 7.051 ± 0.032 ± 0.076 


1.208 ± 0.009 ± 0.020 0.997 ± 0.009 ± 0.013 
1.235 ± 0.010 ± 0.020 0.977 ± 0.010 ± 0.013 
1.243 ± 0.010 ± 0.020 0.966 ± 0.010 ± 0.012 


K 


duscb 
dusc 
dus 


0.737 ± 0.012 ± 0.013 0.841 ± 0.010 ± 0.015 
0.692 ± 0.013 ± 0.013 0.836 ± 0.010 ± 0.015 
0.637 ± 0.013 ± 0.012 0.835 ± 0.010 ± 0.015 


1.141 ± 0.023 ± 0.019 0.942 ± 0.019 ± 0.016 
1.208 ± 0.026 ± 0.020 0.956 ± 0.021 ± 0.018 
1.310 ± 0.031 ± 0.022 1.018 ± 0.025 ± 0.018 


P 


duscb 
dusc 
dus 


0.332 ± 0.010 ± 0.004 0.485 ± 0.009 ± 0.008 
0.333 ± 0.010 ± 0.004 0.494 ± 0.009 ± 0.008 
0.343 ± 0.011 ± 0.004 0.489 ± 0.009 ± 0.008 


1.460 ± 0.050 ± 0.031 1.205 ± 0.041 ± 0.025 
1.481 ± 0.053 ± 0.031 1.172 ± 0.043 ± 0.025 
1.427 ± 0.051 ± 0.028 1.109 ± 0.040 ± 0.021 




duscb 
dusc 
dus 


7.077 ± 0.031 ± 0.071 8.573 ± 0.026 ± 0.086 
6.773 ± 0.033 ± 0.068 8.560 ± 0.026 ± 0.086 
6.654 ± 0.034 ± 0.067 8.559 ± 0.026 ± 0.086 


1.211 ±0.006 ±0.014 1 
1.264 ±0.007 ±0.015 1 
1.286 ±0.008 ±0.015 1 


Mercedes events 


Kind 


Quark 


Nquark Nciuon 


Rx R'x 


TT 


duscb 
dusc 
dus 


6.973 ± 0.078 ± 0.085 8.962 ± 0.133 ± 0.096 
6.735 ± 0.084 ± 0.076 9.002 ± 0.134 ± 0.133 
6.700 ± 0.088 ± 0.072 9.028 ± 0.134 ± 0.133 


1.285 ± 0.024 ± 0.022 0.998 ± 0.023 ± 0.012 
1.337 ± 0.026 ± 0.029 0.988 ± 0.023 ± 0.012 
1.347 ± 0.027 ± 0.029 0.974 ± 0.024 ± 0.012 


K 


duscb 
dusc 
dus 


0.862 ± 0.025 ± 0.014 0.978 ± 0.041 ± 0.016 
0.819 ± 0.027 ± 0.013 0.982 ± 0.042 ± 0.011 
0.773 ± 0.027 ± 0.013 0.981 ± 0.042 ±0.011 


1.135 ± 0.058 ± 0.021 0.881 ± 0.046 ± 0.018 
1.199 ± 0.064 ± 0.017 0.886 ± 0.049 ± 0.016 
1.268 ± 0.070 ± 0.018 0.917± 0.052 ± 0.017 


P 


duscb 
dusc 
dus 


0.401 ± 0.022 ± 0.006 0.656 ± 0.040 ± 0.013 
0.422 ± 0.023 ± 0.012 0.636 ± 0.038 ± 0.013 
0.408 ± 0.025 ± 0.013 0.674 ± 0.041 ± 0.016 


1.635 ± 0.134 ± 0.049 1.269 ± 0.106 ± 0.046 
1.507 ± 0.123 ± 0.098 1.114 ± 0.092 ± 0.065 
1.650 ± 0.142 ± 0.115 1.192 ± 0.104 ± 0.071 




duscb 
dusc 
dus 


8.467 ± 0.066 ± 0.098 10.91 ± 0.113 ± 0.110 
8.085 ± 0.071 ± 0.097 10.94 ± 0.113 ± 0.123 
7.942 ± 0.074 ± 0.089 10.99 ± 0.114 ± 0.124 


1.288 ±0.017 ±0.018 1 
1.353 ±0.018 ±0.025 1 
1.384 ±0.019 ±0.024 1 



Table 3: Multiplicities and ratios of multiplicities of identified particles in quark and 
gluon jets . 

are produced in preference. All positively charged protons were omitted in the 
corresponding momentum range {p < 2.7 GeV/c) to study this effect. 

3. Particle identification 

To take uncertainties of the particle identification into account, the results for dif- 
ferent particle identification cuts (loose, standard, and tight) were compared. 

4. Purity correction of the jets 

The flavour composition of the normal mixture sample has been varied by imposing 
cuts of different strength to the event sample. In this way three samples were 
obtained, one with the flavour mix of Z decays, one which was depleted in b events 
and one depleted in b and c events. The first and second sample were used to obtain 
pure "duscb" and "dusc" results and the third sample to obtain "dusc" and "dus" 
results using the Monte Carlo. 

Furthermore the results were compared to those obtained by using the Cambridge 
algorithm |^ instead of the Durham algorithm. The change of the multiplicities then 
is typically 2%; changes of the ratios and double ratios are much smaller. Finally a 
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Data 


JT 74 def. bary. 


JT 74 


JT 73 


AR 


HW 


Y Events 


K+ 


0.997 ±0.009 ±0.013 


0.98 


0.99 


1.01 


1.00 


1.05 




0.942 ±0.019 ±0.016 


0.94 


0.89 


0.82 


0.88 


0.69 




1.205 ±0.041 ±0.025 


1.29 


1.58 


1.41 


1.49 


0.94 


Mercedes Events 


K+ 


0.998 ±0.023 ±0.012 


1.00 


1.01 


1.02 


1.01 


1.05 


R'k+ 


0.881 ±0.046 ±0.018 


0.94 


0.90 


0.81 


0.88 


0.70 


R'p 


1.269 ±0.106 ±0.046 


1.20 


1.43 


1.38 


1.37 


1.11 



Table 4: Normalized multiplicity ratios R'x (for dusch quarks) compared to the predic- 
tions from the Monte Carlo simulations (JT 74 def. bary. = Jetset 7.4 PS with default 
baryon production, JT 74 = Jetset 7.4 PS, JT 73 = Jetset 7.3 PS, AR = Ariadne 
4.08, HW = Herwig 5.8C). 



Particle 


This Paper Y 


Delphi 18 


OPAL three-jet 14 


Argus @ 


pp 


0.997 ±0.009 ±0.013 
0.942 ±0.019 ±0.016 
1.205 ±0.041 ±0.025 


0.930 ±0.040 ±0.020 
1.120 ±0.110 ±0.040 


1.016 ±0.010 ±0.010 
0.948 ±0.017 ±0.028 
1.100 ±0.024 ±0.027 


1 (def.) 
0.86 ±0.31 
1.58 ±0.10 



Table 5: i?V from measurements of different collaborations. 



systematic error of < 2% due to track reconstruction losses as determined from the 
overall multiplicity measurements is assumed. As both systematic errors discussed 
apply to particles in general, they are expected to cancel in the ratios and are therefore 
not included in Table ^. 

All systematic errors discussed above apply to spectra of identified particles. However, 
the statistical uncertainty here in is general much bigger than the systematic error due 
to the binning of the data. Moreover many systematic uncertainties will cancel in the 
gluon to quark ratios which are the main subject of this paper. Systematic errors have 
therefore been neglected in the errors shown in the particle distributions. 

3.2 Momentum Spectra 

Figure ^ shows the momentum spectra of identified hadrons in quark {dusch) and 
gluon jets for Y events. The momentum spectra of kaons and protons differ significantly 
from those of pions. Pions are produced mainly at low momentum, both in quark and 
gluon jets. The likely explanation is that pions are often low energy decay products 
of unstable particles. The Monte Carlo generators Jetset, Ariadne, and Herwig 
describe the gross features of the measured Delphi data. The momentum distribution 
of kaons in gluon jets is best described by Ariadne. The Herwig model shows a 
considerable weakness concerning the description of kaon momentum spectrum in gluon 
jets. The multiplicity of fast kaons is clearly underestimated. The momentum distribution 
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1.11 
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0.36 


0.99 
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Tabic 6: Systematic errors for A^^, A^^, Rx and B!^ in Y events. Here X denotes the 
particle species tt, p or all charged particles. The summed error displays the quadratic 
sum of the individual errors. 
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Figure 5: Momentum spectra of identified hadrons in quark and 
gluon jets a)-c) spectra of pions, kaons, and protons in quark 
jets; d)-f) corresponding spectra for gluon jets in events with 
Y topology. The predictions of the generator models Jetset, 
Ariadne und Herwig are drawn as hues. 



Gluon/duscb (g/duscb)p / (g/duscb), 



+ 

n 


• DELPHI + 

— Jetset 7.4 

— Jetset 7.4, def. bary. prod. 
••■ Herwig 5.8C 


^) 




• e) 




■••■f) 



1 10 1 10 

p [GeV/c] 



Figure 6: Ratios of the momentum spectra of identified hadrons 
in gluon and quark jets of Y events; a)-c) ratios of the spectra 
of pions, kaons, and protons in gluon jets to those in quark jets; 
d)-f) corresponding spectra normalized to the ratio gluon/quark 
for all charged particles. The predictions of the generator mod- 
els Jetset, Jetset with default baryon production model and 
Herwig are drawn as lines. 
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of protons in gluon jets is well modelled by the Jetset and Ariadne generators but not 
by the Herwig model. 

Figure ^ shows the ratios of the momentum spectra of identified hadrons in gluon 
and quark jets . This measurement is an improvement of our previous publication |T8| . 
More low energy particles are produced in gluon jets than in quark jets for all kinds of 
particles. At high particle momenta this structure is inverted. Figure ^(d,e,f) shows the 
corresponding ratios of the momentum spectra of Figure |^(a,b,c) normalized to the ratio 
of the momentum spectra of all charged particles in gluon and quark jets . 

Figure ^(f) indicates that the proton enhancement in gluon jets is bigger than that 
for all charged particles. The overestimate of the proton production ratio by Jetset or 
Ariadne (not shown) is presumably due to an extra suppression of baryon production at 
the end of the string (i.e. for quark jets, see also Section |3.4|) . A much better description 
is obtained using the default baryon production model without this extra suppression. 

However, no direct conclusions concerning the ratios of the multiplicities can be drawn 
from the normalized ratios as a function of momentum, because the shapes of the momen- 
tum spectra of kaons and protons differ significantly from those of pions which dominate 
the all charged particle sample. 



3.3 Rapidity 

Figure ^ shows the rapidity spectra of identified hadrons in quark and gluon jets and 
Figure || shows the corresponding ratios. For all particles there are in the plateau, i.e. at 
low rj, 1.6-2 times more particles in gluon jets. An excess of particles is expected due to 
the higher colour charge of the gluon. At high t], i.e. in the range of the leading particle 
only few kaons and protons are observed in gluon jets. 



3.4 ^-Spectra 

Figure ^ shows the C,p spectra of identified hadrons in quark and gluon jets . The 
Jetset and Ariadne (not shown) models provide a reasonable description over a wide 
range of the C,p spectrum. The maximum height is different for quark and gluon jets 
indicating different particle rates. The point of intersection of the C,p distributions of 
quark and gluon jets for pions and kaons is approximately the same, ^i"^ ~ 1.73. For 
protons the crossing point between the quark and gluon distributions is shifted to higher 

(s) 

momentum at ~ 0.74. Proton production is enhanced in gluon jets, but preferentially 
at high momenta. This can be seen more clearly in Figure |1^ which shows the normalized 
ratio Rp{C,p). It is observed that this ratio is unity within errors at very small C,p (highest 
momenta) and close to unity also at large (small momenta). A strong deviation from 
unity is, however, visible in the intermediate region (0.7 < < 1.85). 

A surplus of baryon production in gluon jets and the observed kinematical properties 
can be qualitatively understood if baryons are directly produced from coloured partons 
or equivalently from a colour string. In a parton shower colour conservation leads to the 
so-called preconfinement property, that is a local compensation of colour charge in space. 
Alternatively the produced colour charges can always be ordered to form continuous 
chains or strings in space time. These strings appear naturally in the Lund fragmentation 



model |g] and in the progenitor model of Feynman and Field |]I9| . A colour string ends 
at the primary quarks produced in the underlying hard scattering but is spanned over 
the corresponding gluons. Hadron production now can be assumed to proceed via a pair- 
creation of a quark-anti-quark (or diquark-anti-diquark) pair and a corresponding string 




Figure 7: Rapidity spectra of identified hadrons in quark and 
gluon jets compared to Jetset in events with Y topology. 




Figure 8: Ratios of the rapidity spectra of identified hadrons in 
gluon and quark jets from Y events; a)-c) ratios of the spectra 
of pions, kaons, and protons in gluon jets to those in quark jets ; 
d)-f) corresponding spectra normalized to the ratio gluon/quark 
for all charged particles; The predictions of the generator models 
Jetset, Jetset with the default baryon production model, and 
Herwig are drawn as lines. 




Figure 9: spectra of identified hadrons in quark and gluon jets 
compared to Jetset in events with Y topology. 




Figure 10: Normalized ratio R'^ = Rp/vch for different composi- 
tions of the quark event sample (Y events). 
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break-up. A single break-up in the vicinity of a primary quark will produce a leading 
liadron, whereas close to the gluon in the centre of the string at least two breaks are 
needed before a hadron is formed. To produce a baryon a production of a diquark-anti- 
diquark pair is compelling. Now it should be noted that in the centre of a string (i.e. 
in the vicinity of the gluon) more possibilities exist which lead to baryon formation. A 
primary diquark-anti-diquark break up as well as a secondary one following a primary 
quark-anti-quark creation leads to baryon production (see Figure [TT|). 



I Z 
I 




Figure 11: Different possibilities of baryon production in strings. The 
single points denote quarks and the double points diquarks. Open 
points stand for quarks and filled points for anti-quarks. Line a} il- 
lustrates a primary splitting into diquark-anti-diquark in the center of 
the string. Line b} shows the possibilities for secondary diquark pro- 
duction. Line c} shows a diquark-anti-diquark splitting at a string end. 



The latter process (marked b} in Figure |Tl|) may happen in any of the two remaining 
strings similarly to both original endpoints of the string (see Figure ^ c}). The first 
production mechanism (marked a} in Figure [l^) is missing at the string end. This leads 
to the excess of baryons in gluon compared to quark jets. Here it is likely that two 
leading baryons are produced which take a large fraction of the gluon energy. Thus 
it is expected that the excess of baryon production centers at comparably large scaled 
momentum which is indeed observed (see Figure rU). At small momentum, i.e. in the 
momentum range where baryons from the inner part of the string between the jets are 
expected to contribute, the relative portion of baryon production in quark and gluon jets 
is approximately equal. 

Although the above discussion centers around the string model it is based on quite 
general topological properties and is a strong indication that baryon production, and 
presumably also meson production, happens directly from colour objects and an inter- 
mediate step of colour and baryon number neutral objects is avoided. In particular this 
is also indicated by the failure of the HERWIG model to describe the surplus of proton 
production observed in the data. 

In detail the above described mechanism will be complicated by the abundant pro- 



duction of resonant baryons or equivalently the so-called popcorn- mechanism 0. 

Further support to this interpretation comes from the observed strong energy dependence 
of the surplus of baryon production (compare the Argus measurement at ^/s ~ lOGeV 
to this result in Table ^. The surplus of baryons in gluon jets is due to the leading 



18 



baryons. As energy and thus the muhiphcity ratio in gluon to quark jets increases |T[ 
this excess is less and less important in the double ratio R^. The Lund model qualitatively 
describes the decrease of i?^ shown in Table |. 
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Figure 12: Fitted maxima of the spectra. Here k = Ej^t sin y and k = Ebeam denote 
the underlying jet scale. Observe the offset of 0.5 units for the ^* values for the K^. 
The solid lines represent the fit to Y and Mercedes results, the dashed lines represent 
the fit including the e^e~ data. 



Figure |I2| shows position of the maximum, ,^*, of a simple Gaussian fitted to the 



distributions in dependence of the scale n = Ejef sin Omm/I with Omin being the angle with 



respect to the closest jet, here Omin = (for a detailed discussion of jet scales see [0). 
The fit results for ^* are given in Table |^. The maxima of the distributions for protons 
and kaons for quark jets are shifted to smaller values (i.e. higher momenta) compared 
to pions as has been observed previously. It is clearly observed that the maximum is 
shifted to smaller ^p for i^'s in quark compared to gluon jets. To a lesser extent this is 
also observed for vr's and p's. In MLLA/LPHD the maxima of the spectra for gluons 
and quarks are expected to be almost identical ||T0| , p2| . A natural explanation for the 
observed difference especially for i^'s is the leading particle effect. In it has been 
shown that in a scaled momentum range corresponding to the smallest ^-bin in Figure |^ 
leading charged kaons (i.e. those containing a primary produced quark flavour) are about 
three times more often produced than non-leading kaons. 

For all particles the ^* values are bigger for Mercedes than for Y events, i.e. a scale 
evolution of the ^p spectra is observed (see Figure and Table 0). Assuming a general 
linear increase of ^* with the logarithm of the scale k, i.e. 



/(e;)=a + Mn(«) 
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Y events 
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par- 
ticle 




par- 
ticle 




quark 


gluon 


quark 


gluon 


vr 


3.11±0.01 


3.20 ±0.00 




vr 


3.43 ±0.03 


3.49 ±0.02 


K 


2.21 ±0.02 


2.47±0.01 




K 


2.51 ±0.04 


2.77 ±0.05 


P 


2.21 ±0.02 


2.26 ±0.01 




P 


2.52 ±0.05 


2.57 ±0.06 




2.93 ±0.01 


3.05 ±0.00 






3.24 ±0.02 


3.35 ±0.02 



Table 7: Maxima of the distributions. Errors are statistical only. 



the quark jet measurements extrapolate reasonably well to the measurements for overall 
Z events |2^, and for high energy events p5|]. The dashed lines in Figure |12| indicate fits 



including the Z and higher energy data. 

It is further remarkable, that contrary to the predictions of the LPHD model, the 
peak values of the distributions for kaons and protons in quark jets are almost equal 
values for the in Figure |T^ are shifted by 0.5 units to avoid overlaps with the 
proton results.). This observation contradicts the predictions of the LPHD concept that 
the positions of the maxima of the C,p distributions are proportional to the logarithm of 
the mass of the corresponding particle. It has been shown already (see e.g. that 
mesons and baryons show a behaviour which differs from this simple expectation. This is 
a consequence of heavy particle decays and of the partially different masses of the decay 
particles in (predominantly baryon) decays. This statement is qualitatively confirmed by 
this analysis. 



4 Summary and Conclusion 



Based on a sample of about 2.2 million hadronic Z decays collected by the Delphi 
detector at Lep, the production of identified particles in jets initiated by gluons or by 
quarks, was analysed and compared. 

As observed for inclusive charged particles, the production spectrum of identified par- 
ticles was found to be softer in gluon jets compared to quark jets, with a higher total 
multiplicity. The normalized multiplicity ratio {R') for protons in Y events was measured 
to be: 

R'=^ = [y^Z'^'i' = 1-205 ± OMUat ± 0.025,,, 



Np{ch) denotes the number of protons (all charged particles). Herwig underestimates 
both the kaon and the proton production in gluon jets. 

This surplus of baryon production in gluon jets indicates that baryons are produced 
directly from coloured partons or from strings and that an intermediate state of neutral 
clusters (like in the Herwig cluster model) is avoided. This interpretation is supported 
by the scaled energy dependence of the proton excess and by the evolution of the proton 
excess with energy scale. 

Furthermore the ^p and rj distributions were measured and agreement with the Jetset 
and Ariadne models was found. Herwig underestimates both the kaon and the proton 
production in gluon jets. The maxima of the ^ distributions of quark jets, extrapolate 
well with the scale k to those obtained from all events at different centre-of-mass energies. 
For kaons the maximum is shifted to smaller ^* compared to gluon jets presumably 
because of a leading particle effect. 
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